1 

Tradeoff Analysis of Delay-Power-CSIT Quality 



> 

00 



of Dynamic BackPressure Algorithm for Energy 

Efficient OFDM Systems 



(N 

^ ' Vincent K. N. Lau and Chung Ha Koh 

(N 

;_i ■ Dept. of Electrical and Computer Engineering 

< 
(N 

>^ 
o 



Hong Kong University of Science and Technology, Hong Kong 
Email: {eeknlau, eechungha}@ust.hk 



Abstract 

In this paper, we analyze the fundamental power-delay tradeoff in point-to-point OFDM systems under 
imperfect channel state information quality and non-ideal circuit power We consider the dynamic back- 



o 

. pressure (DBP) algorithm, where the transmitter determines the rate and power control actions based on 

the instantaneous channel state information (CSIT) and the queue state information (QSI). We exploit a 
. general fluid queue dynamics using a continuous time dynamic equation. Using the sample-path approach 

- ■ ■ and renewal theory, we decompose the average delay in terms of multiple unfinished works along a sample 

path, and derive an upper bound on the average delay under the DBP power control, which is asymptotically 
accurate at small delay regime. We show that despite imperfect CSIT quality and non-ideal circuit power, 
the average power (P) of the DBP policy scales with delay (D) as P = 0{D cxp{l/ D)) at small delay 
regime. While the impacts of CSIT quality and circuit power appears as the coefficients of the scaling law, 
they may be significant in some operating regimes. 

"This work has been supported by Huawei Technologies. 



I. Introduction 

There is a growing awareness of energy efficiency of the wireless infrastructure. In |IT], ||2], the authors 
considered adaptive power control to optimize an energy efficiency metric, namely the Joule per bit. This 
metric measures the average cost (i.e. energy) to the generated utility (i.e. information bits) and allows for 
accessing the energy efficiency at full loads. However, this metric does not incorporate the delay aspect. 
Furthermore, in all these works, perfect channel state information (CSIT) is assumed and the issues of packet 
errors are ignored. There have been several works that consider robust power control by taking into account 
of the imperfect CSIT. In ||3l, ||4l, the authors derived a power adaptation algorithm to maximize the system 
goodput accounting for potential packet errors due to imperfect CSIT. There are also some works that consider 
MIMO design with imperfect CSIT |[3l, Q or limited feedback |[6l. While these works have dealt with the 
impact of imperfect CSIT, they have ignored the burstiness as well as the delay performance of the source 
data. 

In general, it is quite challenging to consider delay optimization in wireless systems because that involves 

a joint consideration of both the information theory (to model the PHY dynamics) and the queueing theory 

(to model the delay dynamics). One general approach to delay optimal control is to use Markov Decision 

Process (MDP) and the optimal control is given by the Bellman equation iQ. However, it is well known that 

there is no simple solution to the Bellman equation and one cannot obtain viable solutions using brute force 

value iteration or policy iteration. Another approach is to adopt Lyapunov theory to derive the throughput 

optimal control policju ||9l. While the control policy derived is adaptive to both the CSIT and QSI, the delay 

performance of such schemes are not fully understood. In flQl, the authors derived the dynamic backpressure 

(DBP) control and the well-known asymptotic power-delay tradeoff of 0{l/V), 0{V) as F — )• oo, i.e., at 

large delay regimeB Furthermore, all these existing works have assumed perfect CSIT and it is not clear how 
'Throughput-optimal policies are the set of policies which can stabilize the queues in the system if the arrival rate vector is inside 
the stability region (Sj. 

is a constant which determines the tradeoff between power consumption and delay performance. We use the following notation 
in 111! to characterize the asymptotic behavior a function g{x) as x — >■ x*: g{x) = 0{f{x)), if lim sup^^^. < oo; g{x) = 
r!(/(a;)), if limsup^^^.^M < oo; g{x) = Q{f{x)), if g{x) = 0{f{x)), and g{x) = n{f{x)). 



the CSIT quality will affect the underlying delay power tradeoff. 

In this paper, we are interested to study the inter-relationship between energy efficiency, CSIT quality and 
delay performance of the DBP algorithm for OFDM systems. The following are some first order technical 
challenges that have to be overcome: 

• Coupling between the Control Policy and the Queue Dynamics: The DBP algorithm is adaptive 
to both the QSI and CSIT and hence, this introduces coupling to the queue dynamics and the control 
actions at each frame. This coupling makes the delay analysis of the system extremely difficult because 
there is no closed-form expression for the steady state distribution of the queue length for such dynamic 
policy. While there are some works in the literature that analyze the delay-power tradeoff for DBP 
algorithms using Lyapunov bounds ifTOl . such analyses are focused on the first order analysis. However, 
such Lyapunov bounding techniques are usually very loose and cannot be used to study the second order 
impacts of system parameters such as how the CSIT quality affects the delay performance. 

• Coupling of Imperfect CSIT and Energy Efficiency: For a given delay requirement, we should 
reduce the transmission data rate so as to increase the transmission time ifTlll . |[T3l for energy efficient 
communication^ However, the underlying tradeoff changes when circuit power and imperfect CSIT 
are taken into account. By circuit power, we refer to the power consumption in the RF transmission 
chain, which can be assumed to be constant irrespective of the transmission rate. Thus, it is not always 
energy efficient to lengthen the transmission time up to the delay limit because the energy expenditures 
resulting from circuit power are proportional to the transmission time. On the other hand, due to imperfect 
CSIT, the transmission power required to support a certain goodput depends on the CSIT quality. Yet, 
quantifying this relationship for OFDM systems with joint encoding across subcarriers is not trivial 
because that involves finding the CDF of the mutual information. 

• Performance in Small Delay Regime: Another potential limitation of the Lyapunov bounding technique 
is that it works for asymptotically large delay regime, which is usually not the regime we are interested 

in. There are not many works on the analysis of small delay regions. In ||TT|. the authors studied the 

^This is because the transmission power increases exponentially with data rate. 



delay-power tradeoff of a bounded rate scheme in fading channels at small delay regime under perfect 
CSIT. However, the delay performance of DBP algorithm in OFDM systems at small delay regime as 
well as the impact of CSIT quality are still not well-understood. 
In this paper, we overcome the above challenges using continuous time stochastic calculus |[T4l . |[T5l . 
For instance, we introduce a virtual continuous time system (VCTS) and model the fluid queue dynamics 
using a continuous time dynamic equation |[T3l . |[T6l . Using a calculus approach, we first derive closed-form 
expressions for the unfinished work in between arrivals. To analyze the end-to-end average delay performance, 
we adopt and extend the sample -path approach ||T5l and decompose the average delay in terms of multiple 
unfinished works along a sample path using renewal reward theory. As such, we obtained an asymptotically 
accurate upper bound on the average delay under the DBP power control at small delay regime. We show 
that despite imperfect CSIT quality and non-ideal circuit power, the average power (P) of the DBP policy 
scales with delay {D) as P = 0{D exp{l/D)) at small delay regime. While the impacts of CSIT quality 
and circuit power appears as the coefficients of the scaUng law, they may be significant in some operating 
regimes. 



TABLE I 

Notation used throughout the paper 



Symbol 


Meaning 


At 


time duration of a scheduling slot 


T 


time duration of a arrival period 


B 


the average of aixiving packet size 


Uik) 


total number of remaining bits at fc-th scheduling slot in a queue 


U{t) 


fictitious queue state at continuous time t 


£ 


target packet error rate 


np 


the number of subcamers 




the variance of the CSIT error 




per-period average unfinished work 




per-period average energy consumption 


ju 


continuous time per-period average unfinished work 




continuous time per-period average energy consumption 


d(n) 


average (end-to-end) delay of a policy Q 


m) 


average power consumption of a policy 



II. System Model 

In this section, we shall elaborate the system model of the OFDM link, including the physical layer model, 
the bursty source model, the queueing dynamics and the power consumption model. 

A. Frequency Selective Fading Channel Model and the Imperfect CSIT Model 

We consider a point-to-point OFDM system with np subcarriers. The number of resolvable paths in the 
frequency selective channel is given by = -^j- , where W is the signal bandwidth and A/c is the 
coherence bandwidth. The channel impulse response can be described by: 



h{t,v) 



where ^ is the time delay of the l-th path and hi{t) is the corresponding circularly symmetric complex 
Gaussian (CSCG) random fading coefficients with zero mean and variance af (af defines the power-delay 
profile). Using niT'-point IFFT and FFT in the OFDM system, the received signal in the frequency domain 
is given by: 

where Sn and Zn are the transmit and receive signals, respectively, of the n-th subcarrier and Wn is the 
i.i.d. complex Gaussian noise with zero mean and normalized variance 1/np (so that the total noise power 
across the np subcarriers is unity). Note that Hn = J^iJi^ie '^p for all n, which is the FFT of the 
time-domain channel fading coefficients {ho, - ■ ■ ,hN^-i}. 

For simplicity, we consider a TDD system and the transmitter obtains an estimate of the CSIT based on the 
reciprocal reverse channel IITtI . However, due to the channel estimation noise as well as the TDD duplexing 
delay, the estimated CSIT may be outdated. Assume that the CSIT is estimated using MMSE prediction in 
the time domain, the CSIT model in the time domain is given by: 

hi = hi + Ahi, Ahi^CAf{0,ali), I G {0,1, ...,Nd-l}. 

where cj^ ^ = 1 — -g^?^ Jo(27r//)T), Ep is the pilot SNR, Jq is a Bessel function of the first kind of order 0, 
fi) is the Doppler shift, and r is the duplexing delay lITSl . |[T9ll . Thus, the estimated CSIT in the frequency 



domain Hn after ni^-point FFT of {hi , ...,hNa~i} is as follows: 

Hn = Hn + AHn (1) 

where Hn is the actual channel state information (CSI) of the n-th subcarrier and AHn represents the CSIT 
error. The CSIT errors Ai7„ is CSCG with zero mean and variance = i' correlation of the 

CSIT error between the ni-th and n2-th subcamers is given by: E[Ai?„j Aiif^^ ] — '^1=1^ h 

We model the packet error solely by the probabiUty that the scheduled data rate exceeds the instantaneous 
mutual information. Note that the packet errors due to imperfect CSIT is systematic and cannot be eliminated 
by simply using strong channel coding. Therefore we shall exploit diversity to protect the information 
from channel outage to enhance the chance of successful delivery to the receiver. Specially, the encoded 
symbols are transmitted over the frequency domain via a random frequency interleaver. The conditional 
packet error probability (PER) e (conditioned on the CSIT H and the QSI U) of a transmission with data 



rate r (nat/sec) is given by: Pr 



H , U 



where Ptx is the transmit power. 



H = (^Hq , ... , Hrip-i^ indicates the estimated CSIT which is described in (dJ and U denotes the queue 
length (in number of bits) {U is defined in Section II.B]. 

Remark 1 (Power Allocation over Subcarriers): In this paper, we have assumed uniform power allocation 
over subcarriers for the following reasons. First, there is no known closed form expression conditional 
cdf of mutual information X]n=o^ + ^'^1^"^ )• Second, for moderate CSIT quality, the performance 

bottleneck is the packet errors (due to the residual uncertainty of mutual information given inaccurate CSIT). 
As a result the first order factor is diversity. Based on uniform power allocation, we have shown (Lemma 
[T]) that full diversity order can be captured. Third, we have compared the performance of optimized power 
allocation (obtained by stochastic gradient method EDI . 11211 ) versus uniform power allocation in Fig. |6]and 
it is shown that power optimization over subcarriers only shows marginal gain at moderate CSIT quality. ■ 



B. Bursty Source Model, Queue Dynamics and Dynamic Rate Control Policy 

In this paper, we consider a bursty bit flow model. The time dimension is partitioned into scheduling slots 
of duration At and indexed by k. Multiple slots are grouped as a frame of duration T and indexed by m as 
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Bi 



Scheduling slot 
duration, At A 



Bm-l 



One period, T 



Fig. 1. The random arrival and delay deadline model. 

illustrated in Fig. [T] We have the following assumptions regarding the CSI and the bursty source model. 

Assumption 1 (Quasi-Static CSI): For notation convenience, we denote H(/c) = i^HQ{k) , ... , Hnp-i{k))^ 
andH(A;) = {Ho{k) , ... ,Hn,-i{k)y as the rii? X 1 dimension CSIT and CSI vectors, respectively, at the 
/c-th scheduling slot. The CSI H(A;) is assumed to be quasi-static within a scheduling slot and i.i.d. between 
scheduling slots. ■ 

Assumption 2 (Bursty Source Model): Let Bm be the random new arrivals (in bits) at the m-th frame. The 
arrival process {Bm} is i.i.d. over m according to a general distribution Pr(S) with average E[;Bm] = TB 
where B is the average arrival rate per second. ■ 

The QSI U (k) is defined as the unfinished work (i.e. the total number of remaining bits) in the queue 
at the beginning of the A:-th scheduling slot. Let x(^) = (U{k), H(A;)) be the system state at the /c-th 
scheduling slot. Given an observed system state x, the transmitter adjusts the transmit data rate according to 
a stationary rate control policy defined below. 

Definition 1 (Stationary Rate Control Policy): Let r{k) be the rate allocation action of the OFDM trans- 
mitter at the A;-th scheduling slot. A stationary rate control policy ^2 is a mapping from the system state x 
to a rate control action r. Specifically, r{k) = for all k. ■ 

Given a stationary rate control policy J7, the queue dynamics is given by: 



U{k + l) = [U{k) -r(A;)(l-e(A:))At]++ S^^j l(mod(A;, T/At) = O) 
where x'^ = max{x, 0} and e{k) G {0, 1} is the packet error indicator at the k-t\\ scheduling sloe 



(2) 



We assume there is an error-free and delay-free ACK/NAK feedback from the receiver to the transmitter 



C. Power Consumption Model 

At the transmitter, the power consumption is contributed by the transmission power of the power ampUfier 
and the circuit power of the RF chains (such as the mixers, synthesizers, phase-lock loop and digital-to-analog 
converters). The transmission power Ptx in general depends on the transmitted data rate r as well as the 
CSIT quality. Specifically, the transmitted data rate is given by r = nirlog(l + ^/(e,(Te,H)) for a target 
PER e |[22ll . where /(e, cTg,H) is a "black box function" which characterize the behavior of the underlying 
PHY under imperfect CSIT. In other words, for a given data rate r, CSIT error Ug and target PER e, the 
minimum required transmission power is given by 



PtJr; H) = ^ 4 (3) 

For uniform power-delay profile^ the following lemma summarizes an asymptotically accurate relationship 
at high and low SNR. 

Lemma 1 (Relationship between Transmit Power and CSIT Quality): Under uniform power-delay profile, 
for a given data rate r, CSIT error al and target PER e, /(e,(Tg,H) in dS]) is given by: 

f{e,alYl) = F-l^.{e) (4) 

where = denotes asymptotic equality for high and low SNR and F^}._,2 is the inverse CDF of the non-central 
chi-square random variable ip"^ with non-centrality parameter s^. The chi-square random variable = 
Xlne/fj l^nl^ '^^d degrees of freedom and variance a'^/N^, where Ib is the set of independent 
subcarriers. The non-centraUty parameter is given by s^{Ib) = J2n£iB l-^^^l^ • ' 
Proof: Due to page limitation, please refer to ll22l for the proof. ■ 
The transmission power Ptx{r; H) is a convex increasing function of the data rate r. On the other hand, the 
CSIT quality affects the Ptx{r; H) via the variance of the inverse chi-square CDF ^^2^^2(2;). Specifically, a 

larger transmission power is required for the same data rate r when the CSIT error increases. Fig. |2]illustrates 

^Uniform power delay profile is known to be the worst case profile in frequency selective fading channels 1231 . 1241 . As a result, 
the closed form expression for /(e,crf,H) in Lemma [T] represents the worst case profile. For general power delay profile, there is 
no closed form expression for /(e, a^, H) but it can be obtained via offline PHY level simulation. 
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3.5 

Rate(nat/sec) 



4.5 



5.5 



Fig. 2. Required transmission power (dB) versus data rate with different CSIT error a1, target PER=0.01, single subcarrier, bandwidth 
= IHz. 

the Ptx{r] H) versus r at different CSIT errors under uniform power-delay profile. Observe that a larger 
transmission power is required for the same data rate r when the CSIT error increases. 

On the other hand, the circuit power Pcct is constant regardless of the transmission data rate. The trans- 
mission circuit is assumed to be on during the transmission of a burst. Hence, the total power consumption 
at the transmitter during a burst is given by: 



Note that there are two conflicting dynamics due to (a) bursty arrivals and (b) energy efficiency regarding 
whether one should increase the transmission time to finish the data backlog. It is known |[25l that one 
should try to reduce the data rate when the buffer level is small so as to decrease the chance of the buffer 
being empty. Furthermore, the transmission power increases exponentially with the data rate and hence, one 
should increase the transmission time so as to save energy |[T3l . On the other hand, when the circuit power 
is taken into consideration, the tradeoff dynamics will be changed because the energy expenditures resulting 
from circuit power are proportional to the transmission time of the system [26]. As a result, it is not always 
advantageous to increase the transmission time. 
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III. Dynamic Backpressure Power Control 

In this section, we first focus on a queue stabilization problem and derive the throughput optimal policy 
DBP using the Lyapunov function L{U) = /2 

A. Preliminaries of Stochastic Stability 

We first introduce a few definitions. We say that the OFDM link with bursty arrival is strongly stable if: 

lim sup — E\U{k)] < oo . 

Definition 2 (Stability Region): The stability region Aq of policy i7 is the set of average arrival rates A 
for which the system is stable under Q. The stability region of the system A is the closure of the set of all 
average arrival rates A for which a stabilizing control policy exists. Mathematically, we have A = U Aq, 
where G denotes policy space. ■ 

Definition 3 (Throughput-Optimal Policy): A throughput-optimal policy dominate^j any other policy in G, 
i.e. it has a stability region that is the superset of the stability region of any other policy in G. Therefore, it 
should have a stability region equal to A. ■ 

In other words, throughput-optimal policies ensure that the queueing system is stable as long as the vector 
arrival rate is within the system stability region A. Note that the throughput optimal policy is not unique. 
While throughput optimality does not guarantee delay optimality in the system, the former policy can still 
improve the delay performance. Furthermore, using Lyapunov analysis techniques 0, the throughput optimal 
policy derived usually has a simple form, which is desirable for implementation. 

Consider the Lyapunov function L{U) = ll27l . and define the one-step Lypunov drift AL{U) as: 

A {L{U{k))) = E [L{U{k + 1)) - L{U{k)) \ U{k) ] . (5) 
The following Lemma summarizes the results on the Lyapunov drift. 

Lemma 2 (Lyapunov Drift): Let A{k) be the arrival process of queue dynamics Furthermore, let 

A{k) < Ajnax and r{k){l — e{k))At < Rmux, for some positive A^ax and Rmax- The one step Lyapunov 
'a policy 0,1 dominates another policy ^2 if Aoj C An-^. 
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drift for the OFDM link with imperfect CSIT is given by: 

AiL{U{m < ^{Al,^ + RlJ-E[U{k){r{k){l-eik))At-A{k)}\Uik)] (6) 
Proof: The proof follows similar technique as in ||9] and is omitted. ■ 

B. Dynamic BackPressure Algorithm 

Based on the Lyapunov drift A {L{U {k))) in Lemma [2l the DBP algorithm can be derived by maximizing 
the negative drift term in Given an observed state x any scheduling time slot, the instantaneous data 
rate r£)Bp(C/, H) is given by: 

rDBP (U, h) = argmax {?7r(l - e)At - V [PUr; H) + P,,t) At} (7) 

where F is a constant which determines the tradeoff between power consumption and delay performance. 
Solving the problem in ^ for the DBP, the instantaneous data rate ruspix) is given by: 



rDBP ix) = np 



log 



(U{l-e)f{e,ai,li) 



+ 



\ V 

During a burst transmission, the instantaneous power consumption under DBP is given by: 



(8) 



U{l-e)nF _ Up I p -f U(l-e)nF _ np ^ r, 

, otherwise. 
Remark 2 (Multilevel Water-Filling Structure of the DBP): The power control action in (|9]l is a function 

of both CSIT and QSI (where it depends on the CSIT indirectly via the noncentrality parameter s^{Ib) = 

J2neiB I "1^ /(^j'^ejH)). It has the form of multilevel water-filling structure where the transmission 

power is allocated according to the CSIT but the waterlevel is adaptive to the QSI. The parameter V acts like 

the Lagrange Multiplier which determines the tradeoff between power consumption and delay. Furthermore, 

the Pcct affects the power control (or rate control solution) in (0 and ([S]) by introducing a penalty proportional 

to a burst transmission time. ■ 

IV. Delay-Power Tradeoff of DBP with imperfect CSIT 

While the DBP is throughput optimal, we are interested in studying the end-to-end delay performance and 
the relationship between the average delay, average power and the CSIT quality. In this section, we shall 
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analyze the power-delay tradeoff using continuous time approximation and renewal process theory. As we 
shall illustrate, this approach not only yields first order tradeoff relationship (at small delay regime) but also 
yields the second-order impacts due to imperfect CSIT and static circuit power Pcct- 



A. Per-Period Unfinished Works 

To analyze the average delay, we first focus on the analysis of one arrival period T. Specifically, define 
the per-period average unfinished work {Uq) and the per-period average energy consumption j^{Uo) as: 



jU{Uo) = E 



A:=0 



U,At U< 



E 



(10) 
(11) 



where Uq is the leftover bits at the buffer at the starting epoch of a period T. To compute j^{Uo) and j^{Uo), 
we shall adopt a continuous time approach. Specifically, we define a virtual continuous time system (VCTS) 
as follows: 

Definition 4 (Virtual Continuous Time Systems): A virtual continuous time system is a fictitious system 
with a continuous queue state U{t) and channel state H(t). The fictitious queue state evolves according to 
the following dynamic equation: 

^ = -E [r* {U{t), H(t)) \U (t)] (1 - e) (12) 

where r* is the DBP rate policy in (HJ and the fictitious channel state is a white procesqj with identical 
distribution as the actual CSIT H. Using the VCTS, the corresponding continuous time average per-period 
unfinished work J^(U,t) and average per-period energy consumption J9{U,t) are defined as follows: 



jU{U,t)=E 

[ /-T 

J9{U,t) =E 
where x{t) = {U{t), H(t)). 



U{s) ds U (t) 



L-Jt 



<7 r* [/(s), H(s) ,xis)) ds U{t) 



(13) 
(14) 



^It means that for any given t, the distribution of random variable H(t) is the same as that of the actual CSIT H. Furthermore, 
we have E[H(t)] = E[H] = 0, E[H(ti)H(t2)] = E[HH^](5(ti - ti). 
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Discrete Time Queue Trajectory 



L'(Cm-l)r)L. 



U{mT) 



.1^ 



( f{U{mT)) 



Continuous Time Queue Trajectory 



U(mT) ' 



Time 




(?JJ-l)r rar-if'' '-mT+At (m+l)r (m + 2)r mT 

U{(m-l)T)^U((m-l)T -&t) + B,_, U(mT)^U{mT - At) + U{(m + r)T)^U {(m + l)r - At) + 



(m + l)r 



Fig. 3. Illustration of the discrete time and VCTS queue trajectory and the unfinished works {j^ , J'^} in one inter-aiTival interval 
T. 

Note that ([T3]).(fT4l) are the continuous-time counterparts of the discrete-time versions in (([TO]l).( (ITTI l). They 
measure the total queue length (total area) of the queue trajectory and total energy consumption (total area) 
of the power trajectory during an inter-arrival interval, respectively. The queue trajectory and {j^, J^} are 
illustrated in Fig. |3] To derive the actual per-period average unfinished work and energy consumption 
and j9 in the discrete time, we first determine the continuous time counterparts using the following lemma. 

Lemma 3 (Verification Lemma of and J9 in VCTS): If there exists a continuous and differentiable func- 
tion iU,t) and J9*(U,t) satisfying the following partial differential equations (PDEs), respectively: 



(1 -e)E r* (C/(t), H(t)) Uit) +U{t)+^ = 



E 



dU 

g{r*i-),x{t)) m 



dJ9 

— 

dU 



dt 



1 -e)E r* {U{t), H(t) U{t) 



dJ9 



(15) 
(16) 



then, {U,t) and J9*{U,t) are the total average unfinished work and total average energy consumption 
of the VCTS under the DBP in (O. ■ 
Proof: Note that (fT5l).(fT6l) resembles the Bellman equation of the discrete time dynamics in ([TOl i. (ITTI ). 
The proof is obtained using Taylor expansion of the value function and using the divide-and-conquer principle 
from ([T3]).(fT4]i. Please refer to Appendix A for details. Using Lemma [3l the areas of the queue trajectory 
and the power trajectory can be obtained by solving the PDEs in (fT5l).([T6l). ■ 
Finally, j^{Uo) and j9{Uo) are related to the continuous time counterparts by the following Theorem. 
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Theorem 1 (Relationship between the Continuous Time and Discrete Time Unfinished Works): Let {U,t) 
and J9* iU , t) be the solutions of the PDEs ( fTSl ) and ( fT6l ). respectively, in the virtual continuous time system. 
For sufficiently small At, the discrete-time per-period average unfinished work and energy consumption 
and j9 are given by: i^(C/o) = 'W\Uo,Q) + 0(At) and j9(C/o) = J9*([/o,0) +0{M). ■ 

Please refer to Appendix B for the proof. As a result of Theorem [T] we can focus on the continuous time 
equations to solve the area of queue and power trajectories and we can be assured that the solutions obtained 
will be accurate up to 0(At). Using Lemma [3] and solving the associated PDEs in (ITSl) and (fT6l ). we shall 
obtain an asymptotically accurate performance bounds of {U ,t) and J9*(U,t) which is summarized 
below. 

Theorem 2 (Performance Bounds of iU,t) and J^* (JJ ,t) in VCTS): The per-period average unfinished 
work (U,0) in ([T3l) and energy consumption J9*(U,0) in (fT4b are given by 



iU,0) < / y{t; (3)dt 



yjt; (3')nFil - e) 
V 



+ : 



Uonpil - e) 



V 



/(e,^T2,H) 



cct 



Uonpil - e) 
V 



(17) 

dt 
(18) 



,/3= E 



log((l-6)/(6,CT2,H)) 



where /3) = exp |^-/3 + Ei^-^) pi (^log(C/o) + /? j - npil - e)e^t 

and (3' = E ( log((l - e) /(e, o-^, H)))+ . Ei(x) = J^^e^tdt for x > 0, is the exponential integral 
function. The bounds are asymptotically accurate as V approaches 0. ■ 
Please refer to Appendix C for the proof. 



B. Delay Analysis for Deterministic Arrivals 

In this subsection, we establish the relationship between multiple per-period unfinished works and average 
end-to-end delay under the assumption of deterministic arrivals. Specifically, we assume the bit arrival Bm 
is deterministic (given by B). The average bit arrival rate (bits per seconds) is given by B = B/T. Such 
an arrival model embraces VoIP as well as other delay-sensitive source models derived from constant bit 
rate (CBR) encoders |[28l . l|29l . Let J(f]*) and g{^l*) be the average end-to-end delay and the average 
power consumption, respectively, under DBP poUcy Q* in actual discrete time systems. d{Q*) and g{0.*) are 
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represented by the combination of per-period average unfinished work over multiple periods. We first have 
the following results regarding the steady state leftover bits of the buffer in VCTS due to the accumulated 
arrivals in the previous arrival periods. 

Lemma 4 (Steady State Leftover Bits for the VCTS): Let Lm be the leftover bits of the buffer in VCTS at 
the end of the m-th arrival period. Given Lq < L*, then we have sup,„gjj+ < L*, and L* > satisfies 
the following fixed point equation: 

Furthermore, the fixed point L* exists and is unique. ■ 
Please refer to Appendix D for the proof. Based on Lemma lU the average delay and power consumption 

for deterministic arrival is given by: 

Theorem 3 (Average Delay and Average Power Consumption for Deterministic Arrivals): For sufficiently 

small At, the discrete-time average end-to-end delay and average power consumption g{Q*) under the 

DBP rj* with deterministic arrivals are given by: 

d{n*)< L^JU*{B + L*,0) +0{At) (20) 
B T 

Vm > ^J^*iB,0) + OiAt) (21) 

where B is the number of bits of an arrival packet in each period. ■ 
Proof: Please refer to Appendix E for the proof, where the relationship between the discrete model and 
the VCTS in Theorem [T] is utilized. ■ 
Corollary 1 (Asymptotic Power-Delay Tradeoff of DBP at Small Delay Regime): For sufficiently small Vj^, 
the asymptotic power-delay tradeoff of DBP of the VCTS is given by: 

I I (22) 

log{BE[f{e, al U)]/V) E[/(e, al H)] / 

g{i}*) = n((^ + P,rf^ 1 (23) 

is a parameter that determines the tradeoff between power and delay in the system. For a given data arrival rate, small V 
corresponds to small delay regime. 
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Proof: Please refer to Appendix F for the proof. ■ 
Remark 3 (Interpretation of Results): Note that from Theorem |3] there is an additional O (At) term in (|20l l 
and ((2T]) accounting for the approximation error between the power/delay of the original discrete time system 
and the VCTS. Yet, to simplify discussion, we focus on the first-order comparisons from the power-delay 
tradeoff of the VCTS in Corollary H] 

• Comparison with CSIT-only policy: The power-delay tradeoff result for the VCTS in Corollary [His 
asymptotically accurate as F ^ and this corresponds to small delay regime. For a given CSIT quality 
and Pcct, the conditional average data rate (conditioned on the queue state U {t)) for CSIT-only policy 
is given by: 

\l-e)f{e,aln) 



I^CSIT — E 



rcsiT{m))\U{t) 



E 



npi log 



(24) 



V I \ 

As a result, the power-delay tradeoff of the VCTS for CSIT-only policy ifTTI at small delay regimqf 

is given by g = 0{exp{l/d)). On the other hand, since jog^jy^^ = ^{V), we have the delay and 

power of the VCTS given by d{n*) = O {j^^^), and g{n*) = Q [ v\og\i/v) ) f™"^ Corollary E 

Furthermore, since it is asymptotically accurate for small V (small delay regime), we can conclude that 

the power-delay tradeoff for DBP in the VCTS is given by ^ = 0{dex.p{l/d)). Hence, compared with 

CSIT-only policy, the power consumption of DBP in the VCTS increases slower as delay d tendo to 

0. Furthermore, we could achieve this superior tradeoff performance even with imperfect CSIT quality 

and non-ideal circuit power {Pcct > 0). 

• Effects of CSIT quality: The penalty of CSIT quality is contained in lE[/(e, dg, H)], which appears in 

the coefficients of the tradeoff equations in (l22l ) and (|23] l. For a given target PER e, a larger CSIT error 

corresponds to a smaller E[/(e, Ug, H)]. Fig. |4] illustrates E[F^2^^2(e)] versus CSIT errors o"g for target 

'The delay expression for CSIT-only policy in 111] is derived using a discrete time approach and is given by d = 0{ ^^^^-^ ) + 1. 
On the other hand, the result in this paper is derived using a continuous time approach (VCTS) and the delay is given by d = 
O jag(l/Y^ ^ + 0{At) (for the actual discrete time system in terms of seconds). Hence, they match each other when expressing in 
terms of seconds. 

"'Note that while the delay of the VCTS can go to zero as V ^ 0, the delay of the actual discrete time system cannot go to zero 
and is given by 0(At) when V — )■ 0. This footnote applies to Corollary |2] as well. 
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D.09 D.I 



Fig. 4. Expectation of F^2.^2 {^) versus CSIT error at different target PER. 



PER 10~^, 10~^ and 10~^ under the uniform power-delay profile. 
• Effects of Pact' From (|23] ). the average power consumption has two components, namely the transmission 
power and the circuit power. The term Pcct/{\og{Bf{e,a1,'H.)/V))) corresponds to the circuit power 
consumption, which increases with the burst transmission time in one arrival period. For small V, the 
burst transmission time decreases in the order of 0(l/(log(i?/(e, cjg, H)/y))). ■ 



C. Delay Analysis for Random Arrivals 

In this subsection, we shall extend the analysis to i.i.d. random arrival process where the bit arrivals 
Bm G [0, -Bmax] is generated by a general distribution 'Pt{B). The main results are summarized below. 

Theorem 4 (Performance Bound of Average Delay and Power Consumption for Random Arrivals): For suf- 
ficiently small At, the discrete-time average end-to-end delay dudi^*) and average power consumption 
gudi^*) for DBP O* under i.i.d. arrival process {Bm} are given by: 



{n*)< ^^E[JU*{B + L*^,,,0)] +0(At) 

9udm>inT9\BM+om 



(25) 
(26) 



where I/^ax given by the fixed point of equation ( fT9l ) (with B = Bmax), Bad is the average bit arrival 
rate (bits per seconds) given by Bad = E[j8]/T and the expectation is taken w.r.t. the i.i.d. arrival process 
{Bm}. ■ 
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Proof: Please refer to Appendix G for the proof. 



Corollary 2 (Asymptotic Power-Delay Tradeoff of DBP for Random Arrivals at Small Delay Regime): For 
sufficiently small V, the asymptotic power-delay tradeoff of DBP in the VCTS under random arrivals {Bm} 
is given by: 



In this section, we shall compare our proposed DBP with two reference baselines, namely the CSIT-only 
control policy (baseline 1) and No-CSIT policy (baseline 2). The baseline 1 policy allocates the rate and 
power to optimize the PHY throughput based on CSIT only. The baseline 2 policy always transmit with 
uniform power and fixed rate. In the simulation, we consider both deterministic and random arrivals. The 
OFDM systems has 1024 subcarriers with total bandwidth lOMHz. The scheduling slot duration At is 5msec. 
We simulate 10^ scheduling slot to evaluate the average power and delay for different parameter V. The 
dashed lines that pass through the simulation point for DBP algorithm represents the analytical results in 
(EH) and (I23]l. 

Fig. |5] illustrates the power-delay tradeoff at different CSIT errors a1 = 0.01, 0.05, 0.1. It can be observed 
that DBP simulation results match the performance bounds in Theorem |3] quite closely. In addition, it is 
obvious that DBP has significant gain compared with the CSIT-only policy and No-CSIT policy. It can be 
observed that the simulation points match with the analytical results very well for small V (which corresponds 
to small delay regime). As the CSIT error a1 gets smaller, the power-delay curve has steeper slope, which 
means better tradeoff. The performance gap between DBP and other policies increases at small delay and 
small CSIT error regime. 

Fig. |7] shows that power-delay tradeoff with different circuit power consumption Pcct- It can be observed 
that the effect of Pcct is significant when Pcct is non-negligible from the total power consumption, especially 




(27) 




(28) 



Proof: The above results can be obtained in a similar way as Corollary [J based on Theorem |4] 



V. Simulation Results and Discussions 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Average Delay (sec) 



Fig. 5. Power-delay tradeoff for OFDM link at different CSIT errors for deterministic arrival. CSIT error variance {al = 
0.01,0.05,0.1), traffic loading B = Ik nats/slot, r=100msec, At = 5msec, P^ct = 0%, nF=1024, Nd =16 and target PER =0.01. 
Note that V is a parameter that determines the tradeoff between power and delay. For example, V = 1,2, 4, 6, 10, 15, 25, 40 at the 
marks along the DBP curves. 





DBP with optimized power allocation 



0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4 0.42 
Average Delay (sec) 



0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4 0.42 
Average Delay (sec) 



(a) al = 0.1 



(b) at = 0.15 



Fig. 6. Average Delay versus power consumption at different CSIT errors for deterministic arrival: traffic loading B = Ik nats/slot, 
r=100msec, At = 5msec, P^ct = 0, nF=1024, Nd =16, and target PER =0.01. 



in the large delay regime. Similarly, the DBP has significant gain compared with the CSIT-only policy and 
No-CSIT policy. 

Fig. [8] illustrates the power consumption versus CSIT errors at different delay requirements. It can be seen 
that when the CSIT error a1 increases, the minimum required power for satisfying the delay requirement 
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Fig. 7. Power-delay tradeoff for OFDM link at different Pcct for deterministic arrival. Traffic loading B = 0.8k nats/sec, T=100msec, 
At = 5msec, Pcct = [0,7,10]dB, ac = 0.05, nF=1024, Nd =16 and target PER =0.01. Note that V is a parameter that determines 
the tradeoff between power and delay. For example, V — 0.5, 1, 2, 4, 6, 10, 15, 25, 40 at the marks along the DBP curves. 
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(a) Deterministic arrivals 



(b) Random arrivals 



Fig. 8. Power consumption versus CSIT errors at different delay requirements: (a) Deterministic arrivals; (b) Random arrivals, traffic 
loading B = 10k nats/slot, r=100msec. At = 5msec, Pcct = 10%, nF=1024, Nd =16, and target PER =0.01. 



increases. In addition, the performance gain of the DBP decreases as the CSIT error increases. 



VI. Conclusion 

In this paper, we consider a tradeoff of power-delay in point-to-point OFDM systems with imperfect CSIT 
and non-ideal circuit power. Using Lyapunov optimization framework, we derive a dynamic backpressure 
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algorithm (DBP), which adapts the rate and power based on the instantaneous CSIT and QSI. To study how 
the CSIT quaUty and circuit power affects the power-delay tradeoff, we introduce a virtual continuous time 
system and derived an asymptotically accurate the power-delay bounds at small delay regime. We show that 
despite imperfect CSIT quality and non-ideal circuit power, the average power of the DBP policy scales with 
delay {D) as 0(1) exp(l/D)). The impact of CSIT quality and circuit power appears in the coefficients of 
the scaling law. 

Appendix A: Proof of Lemma[3] 

In order to obtain the continuous-time area under the queue trajectory in the VCTS, we shall use the 
principle of divide -and-conquer from the definition of J^{U,t) in (fT3] ) and J9{U,t) in ([T4b . Specifically, 
we have 



ju(u,t 



E 



U{s)ds + jU{U,t + tA) U{t) 



Using Taylor expansion J^{U,t + t^) = {U ,t) - - e)r*(-)tA + + 0{tA^). For small 

value of tA, U (s) is assumed to be fixed to U {t) for s e [t,t + 1^]- Thus, removing J^{U, t) on both sides, 
we have 

= E 



U{t)tA - ^(1 - e)r*{-)tA + -^tA + 0{tA') U{t) 



Dividing by Ia and taking the limit tA — gives ([TSl l. (fT6l ) is also obtained by a similar procedure from 
the definition of J^([7, t) in (fT4l ). 



Appendix B: Proof of Theorem [T] 

Similar to Appendix A, using the principle of divide and conquer from the definition of {Uq) in ([TOl l 
and j3{Uo) in ([TT]) . we shall have the following recursive equations in discrete-time systems: 



jU{Uk) = E UkAt + jU{Uk - r*{k){l-e)At)\Uk 



j9{Uk)=E[g{r{k) ,xik))At + j9 {Uk - r*ik)il - e)At) \Uk] , k = 0, 



T 
'At 



(29) 



1 (30) 



where r*{k) is given by the DBP in ([8]). Since {j^, J^} and {J^ , J^} satisfy the discrete time and continuous 
time recursive equations, respectively, we only need to show they are different in 0{At). 
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We then discuss the property of {U,t). Note that {U,t) satisfies the PDE in ([T5l), multiplying At 
in both sides of dTSl ). we have: 



E 



f)jU f)jU 
U{t)At+ ^{-r*{-){l-£)}At+ ^At 



Uit) 



0. 



(31) 



For simpUcity, let r* = E[r*(U{t), H(t))|C/(t)] . Since r* = if U{t) = 0, by Taylor expansion on 
1^0 - r*(l - e)At, t + At), we have: 

Ju0 _ _ e)At,t + At) - J^{U,t) = ^(-^(1 - e))At + ^At + 0{At^). (32) 



By substituting (I32j) into (1311) : 



jU(U,t)= K U{t)At+ JU{U -r*{l-£)At, t + At) U{t) - 0{At^). 
Let t = kAt, U{kAt) = Uk and J^{Uk) = {U , kAt), i.e., sampling at time t = kAt, 



(33) 



jf([/fc) = E UkAt+ Jjl^^iUk+i) 



Compare to the discrete time recursive equation in (|29] . jj^{Uo,0) satisfies it up to 0{NAt'^) = 0{At). 
As a result, the solution of continuous time PDE {Uo,0) will be different from the actual j^{Uo) by at 
most 0{At), i.e. iUo) = (t^0;0) + 0{At). The case of the average energy consumption is obtained 
similarly. 

Appendix C: Proof of Theorem [2] 

In the proof, we shall fist derive the bounds for the average departure rate and power consumption of 
the DBP policy, which is asymptotically tight at small delay regime (small V). Based on these bounds, we 
can derive the upper and lower bound of the queue trajectory and the corresponding bound of the average 
unfinished work and energy consumption. 

First of all, we have the following lemma on the conditional average rate and power of the DBP policy. 

Lemma 5 (Bounds on Conditional Average Policy): The conditional average rate and transmission power 
of DBP can be bounded by: 



fDBp{U{t)) = E \rDBp{x)\U{t)\ £ '^low 

{Uit)), rup{U{t)) 



(34) 
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Queue length U 



Queue length U 



Fig. 9. Comparison of the actual and approximated rate and power control policies. 



9DBpiUit))=K gDBpix)\U{t) > gi^jUit)) 



(35) 



where 



riUU{t)) = E 



np log 



luit){l-e)f{e,alU) 



rup{U{t)) = n^E 
and^,„^(c7(t)) 



U{t){l-e)n F , _ 
V 



V 



(l-£)/(6,a,^H) 

V 



+ 9dbp\'~^ l*ojj V 



U{t) 



nF\og{U{t)) +nFl3 



\ogU{t)] + log 



U{t) 



np 



\og{U{t)) 



+ npl3', (36) 



As a result, we can approximate fpiBP and gjjBp using riow{U) and giowi^) with asymptotically small 
approximation errors at small V. Fig. |9] illustrates the accuracy of the approximation. 

Proof: The key proof is to find the bound for the average departure rate. Specifically, from dSjl, using 
Jensen's inequality yields following inequality, we have the lower bound: 



roBP U{t) = E 



np log 



j u{t){l-e)f{e,al 



H) 



U{t) 



\ 

Similarly, using the fact that (x + y)+ < (x+ + we have the upper bound: 

(l-£)/(e,a2,H)'^ ^ 



> rion. U{t) . 



(37) 



roBP U{t) = npE 



log 



{u{t)} 



+ log 



V 



Uit) 



< rup{U{t)). 



Now, we shall show the lower bound of goBpi^i^))- Let /(e, cXg, H) = X and q{x) to be the pdf of the 
random variable X, where the randomness is induced by the CSIT H. From we have: 



gDBp{U{t))=¥. 



U(t)np(l — e) np 

V IT 



cct 



Uit) 



I Unp{l — e) np 



V \ 



V 



X 



+ Pcct q{x)dx. 
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9dbp{^{^)) is ^ monotonic increasing function of U{t). Differentiating both sides w.rt. U: 



ddDBpiU) {l-e)nF 



dU 



V 



q{x)dx< ii_£l!^, Vgi,5p(^)>0. 



As a result, we can construct a lower bound of g^^p by the following ODE: 

dgiowiu) _ 



dU 



ii^^, if [/ e 01 



V 



, otherwise 



(38) 
(39) 



where 0i is <U 



9iow{U) > o|. giow{U) always has steeper slope compared with goBpiU) and together 
with the monotonic increasing property of 'SoBpiU), we could establish the lower bound by solving (|38] l 
and (Ell). ■ 
Secondly, based on Lemma |5l we shall derive an upper bound on (U, t) by solving the PDE in Lemma 
[21 Specifically, the queue dynamics U (t) satisfies the ODE in the VCTS in ([T2l) : 



^ = -rDBp{U){l -e)< -riUU)(.l - e). 



(40) 



Hence, the upper bound queue trajectory is given by: 



U{t) < Uup{t) = exp [-/3 + Ei(-i) [Ei i\og{Uo) + - (1 - e)nFeH\\ = y{t; (3), (41) 



Based on ( fT3l ). we have 



J^*iU, 0) «E 



U{t)dt 



< I Uup{t)dt 



which yields the upper bound of ([/, 0) in ([TT] ). 

Finally, we shall derive a lower bound for j9*(C/,0). Using rup{U{t)), we can construct a lower bound 
trajectory by solving = —rup{Uiow)0- ~ ™d the solution is given by: 

U{t) > UioUt) = exp + Ei(-i) [Ei (log([/o) + P') - (1 - e)n,.e^'t]] = /?'), (42) 

At any time i G [0, T], we have 



gioMoUt)) < 9iov,{u{t)) < gDBp{u{t)). 



(43) 
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Hence, 

J9*{U, 0) = / gDBpiU{t)) dt > / 
Jo Jo 



Uiow{t)nF{l - e) ~ Uonpil-e) 

+ 9dbp\Uo) 



dt 



which yields the lower bound of per-period average energy consumption in dTSl ). 

Appendix D: Proof of Lemma H] 

Using the theory developed for the continuous time model in Section IV-A, it is enough to show that Lm 
for the VCTS is bounded by L* for all m. First of all, we shall show that L* is unique. 
Specifically, the leftover bits at the end of the m-th period is: 

^ Lyti ^ lm (44) 

where l„i and are the leftover bits of m-th period of approximated queue trajectory Uup{t) and Uiow{t) 
in (|4TI ) and (l42l ). respectively. Recall that Uiow{t) <U{t)< Uup{t). Since the unfinished work at the start 
epoch of m-th period is the summation of arriving bits B and the leftover bits of the previous (m — l)-th 
period, lm and lm satisfy the foUowings: 

lm = f {B + Im-i) and ljn = f_(^B + Im-l ) 

where 7(x) = exp [-/3 + Ei^-^) [Ei(log(2;) + ^) - npil - e)e^T]] , and /(x) = exp[-/3'+Ei(-i) [Ei(log(x)+ 
/?') — ni?(l — e)e'''T]]. Both f{x) and f{x) are the increasing functions of x since the exponential integral 
Ei is increasing function. Note that the slope of f{x) is given by: 

djjx) _ Ei(-i)[Ei(log(x) + /3)-(l-e)n^e^r] for all x > 0. 

dx log(x) + /3 ' 

Therefore, there exists a unique crossing point between y = x and y = f{B + x) and this proved the existence 

and uniqueness of L* . 

Finally, we shall try to prove that sup^gK+ Lm is bounded by L*. We first claim that lm < L* for all m. 
From ([T9l ) and lm = f {B + lm-i) , we have: 

Ei (log(L*) + ;3) = Ei {\og{B + L*) + /?)-(!- e)nFe^T (45) 

Ei (log(I„) + /3) = Ei {\og{B + Im-i) + /3) - (1 - eWe^T (46) 
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Subtracting (|46] l from (|45] l. we have 



log(L*)+/3 gx /■log(B+L-)+/3 



= / — dx. (47) 

iog(I„o+/3 ^ iiog(B+I,„_i)+/3 a; 

Note that e^/x is positive for x > 0, and log(Im) + /3 > due to /(x) > e~^. If {Im} is not bounded by 
L*, there exists m' such that /m'-i < L* < Im' and S + Im'-i < i? + -L* < -B + Im'- This m' makes the 
RHS of (l47l) positive while the LHS of 1471) becomes negative which means a contradiction. Thus, {Im} is 
bounded by L*. As a result, sup^£]g+ Lm is bounded by L*. 

Appendix E: Proof of Theorem [3] 

In the proof, we shall use Little's law f30l to derive the average delay and power consumption for the real 
discrete time system. 

Specifically, under the stationary DBP policy, the system state x{k) evolves as an ergodic Markov chain 
and hence, there exists a steady state distribution tt^ such as vr^ (xo) = linifc-j^oo Pr [x(^) = Xo] ■ Using 
Little's law and the characteristic of ergodic chain, the average end-to-end delay under the stationary DBP 
policy ri* is given by: 

d{n*) = LE^\U{k)] =L lim y^~^lE[C/(A:)|C/(0)] 
where the E.,^^ is the expectation w.r.t. the steady state distribution of U{k). Using ergodic theory, we have: 

- 1 1 sr^M-l—* 

d{n*)= =- lim — y jU (U^) 

BT M^oo M ^m=0 ^ ^ ' 

where Uq^ is the initial queue length of the m-th period. Note that the unit of i?r is bits, and the unit of 
is bits X seconds, and hence the unit of d is seconds. As shown in Lemma |4l the leftover Lm at the end 
of the m-th arrival period is bounded by L*. Hence, fj™ is upper bounded by B + L*. Using the relationship 
between the unfinished works in discrete time and continues time in Theorem [T] the discrete time average 
delay given by: 

-dm< ±.-jU*(^B + L*)= ^'W\b + L*,0) +0{At). 
Similarly, the average power consumption 'g{^*) is given by: 

g{n*) = l- lim i-y^^"^j9*(C/™) > 1 = 0) +0(At). 
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Note that the unit of T is seconds, and the unit of j9* is Watt x seconds, and hence the unit of 'g is Watt. 

Appendix F: Proof of Corollary [T] 

From the Theorem |3] we only need to obtain (B + L*) and J9*{B, 0) for the continuous time model 
by studying the continuous queue trajectory. 

First of all, we show the asymptotic behavior of L* as V goes to 0. Let x to be /3 for simplifying notation. 
From ([I9]l, we have: L* = e~^+Ei<"'MEi(iog(iJ+L-)+x)-n^(i-e)e1 y g^g^ ^ ^j^^ ^ increase since 

/3 = 0(log(l/y)). For X > 0, increases faster than Ei(x) because 

Ei(x) = logx + 7 + W. < ^.=1 ¥ = ' • 

Hence, as V approaches to 0, we have L* = Q (e~'^) = ©(^j^jy^-^a 

Secondly, we shall obtain the asymptotic area of buffer trajectory. Specifically, we use Fig. [10] to illustrate 
the proof. The upper bound of queue delay can be derived from the area of Uup{t) in (|4T]) . As shown in (l40l) . 
dUup/dt (X log Uup{t), which is decreasing in t. Hence, Uup{t) is a convex function in t and we could upper 
bound J'^ {B + L*, 0) by the summation of the triangle area (A) and the rectangle area (B) as illustrated in 
Fig. [Ma): 

JU*{B + L*,0)< ^{B- La) td + T{L* + La). (48) 
where La = ©(e^'^) and td is the time when Uup{t) = L* + La which is given by: 

td = 7r^4— {Ei MB + L*) + /3) - Ei (log(L* + La) + /?)} = © ^ ^ 



{l-e)nF \log{BE[f{€,alH)]/V), 
From (|48] |. the upper bound of average delay d($7*) is given by (|22]) . 

Then, we focus on deriving J9*{B, 0). Specifically, we use Fig. [TO]to illustrate the proof. We elaborate the 
asymptotic expression using giow{Uiow{t)) where Uiowit) is the lower bound trajectory derived using Tup in 
(l42b . According to ([43] ). ^^^^ is a lower bound of the actual power trajectory and it is a convex function of 
t. As a result, the lower bound of the transmission energy (area of giow) is given by the triangle (C) of the 
tangent line of the 'giowi^) as illustrated in Fig. [TOtb). Let tp to be the time when the tangent line touches 
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Fig. 10. Asymptotic upper bound of the area below the queue and power trajectory. 



zero which is given by 



E 




[ V /(e,a^,H)J 


+ 

+ Pcct 


dt 




U=B 



e 



B 



log(i?E[/(6,a2,H)]/y) 



Note that tp means the transmission time for the burst transmission. On the other hand, we get the same 
order of tp in the case of that "qi^ is a concave function Hence, the asymptotic lower bound for the total 
energy consumption in one period is given by: 



J^\B, 0) > \mon^{B)tp 



np 



V 



~l~ Pcct^p 



which yields the lower bound of average power consumption in 



Appendix G: Proof of TheoremH] 

To derive the bound for the random bits arrival, we shall first derive the upper bound of the steady state 
leftover queue length of a period (where the lower bound is obviously.) 

First of all, from the results of deterministic arrivals in Theorem |3j the average end-to-end delay d{Q*) 
and the average power consumption g{^*) for random arrivals are given by: 



J-E 
BT 



1 



-E 



iB + C,0) + 0{At) 
J^{B + C, 0)1 + 0{At) 



"The tp of the concave case is obtained from gioiu{Uiow{tp)) = 0. 
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where £ is the steady state leftover queue length of a period for random arrivals. From Lemma @] it is 
obvious that L^^^, is the upper bound of £. Then we have the upper bound and lower bound of average 



per-period unfinished work and energy consumption in VCTS: E 



JU*{B + C, 0) < E Ju\b + K,,^,Q) 



and E 



J9*{B + C, 0) > E J9*{B, 0) 



, which yield (1251 ) and (I261 ). respectively. 
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